Anisotropy of ultra-thin ferromagnetic films and the spin reorientation transition 
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The influence of uniaxial anisotropy and the dipole interaction on the direction of the magne- 
tization of ultra-thin ferromagnetic films in the ground-state is studied. The ground-state energy 
can be expressed in terms of anisotropy constants which are calculated in detail as function of the 
system parameters and the film thickness. In particular non-collinear spin arrangements are taken 
into account. Conditions for the appearance of a spin reorientation transition are given and analytic 
results for the width of the canted phase and its shift in applied magnetic fields associated with this 
transition are derived. 



I. INTRODUCTION 

Experimentally it became possible in recent years to 
grow epitaxial thin films of ferromagnetic materials on 
non-magnetic substrates with a very high quality. This 
offers the possibility to stabilize crystallographic struc- 
tures which are not present in nature, and which may 
exhibit new properties of high technological impact. To 
understand the magnetic structure of these systems is a 
challenging problem both experimentally and theoreti- 
cally. 

Generally speaking, for not too thin films the magne- 
tization is in-plane due to the dipole interaction (shape 
anisotropy). However, in very thin films this may change 
due to the increasing importance of surface effects. In- 
deed, at surfaces due to the broken symmetry uniaxial 
anisotropy energies arise which in generally are much 
higher than in the bulk. These anisotropy energies may 
favor a perpendicular orientation of magnetizationtl. Ad- 
ditionally in the inner layers of the film due to strain 
induced distortion bulk anisotropy energies may appear 
absent or very small in the ideal crystal. As a conse- 
quence in these films a reorientation of the spontaneous 
magnetization is observed either as function of film thick- 
ness or as a function of temperature. This spin reorien- 
tatigggguisition has been discussed extensively in the 

Phenomenologically in order to describe the magnetic 
properties, anisotropy coefficients K n compatible with 
the underlying symmetry of the film are introduced which 
are supposed to arise from an expansion of the energy (or 
the free energy at finite temperatures) in terms of the ori- 
entation of the magnetization vector relative to the film. 
These coefficients are then studied experimentally (for 
a review see Ref. ^). In ferromagnetic resonance (FMR) 
experiments, for instance, these coefficients directly enter 
the resonance frequency (for references see for instance 
Ref.|,|). 

Theoretically, it has been shown that the anisotropy 
coefficients K n {T), which are in general temperature de- 
pendent, can be calculated numerically at finite temper- 
atures within mean field theory, starting from a Hamiltpr 
nian with microscopic constant anisotropy parametersQ. 
Furthermore, the temperature dependence of the lowest 
order anisotropy K(T) was determined analytically us- 



ing a combination of-jnean field theory and first order 
perturbation theoryB'EEl. In other approaches the mag- 
netization of the film is calcuJated. directly within mean 
field and spin wave theoryOll3'E3EJ or with-iuil numerical 
calculations like Monte Carlo simulationsEJO. 

In the present paper we describe the ferromagnetic 
film within a classical local spin model with dipolar 
interaction and uniaxial anisotropy. We will concen- 
trate on ground-state properties of thin films in order 
to clarify the discussion and to eliminate all uncertain- 
ties connected with finite temperature calculations. A 
major goal of the present study is the calculation of the 
anisotropy coefficients at zero temperature from the pa- 
rameters of an underlying Hamiltonian. The important 
point is that even in this situation the dependence of 
these coefficients on the microscopic parameters is far 
from being trivial due to non-collinear magnetic states in 
the thin film. It is the purpose of this paper to elucidate 
this behavior. 



II. THE MODEL 

The calculations of the ground-state properties of 
ultra-thin ferromagnetic films are done within the frame- 
work of a classical ferromagnetic Heisenberg model con- 
sisting of L two-dimensional layers with the z-direction 
normal to the film. The Hamiltonian reads 



(ij) 



3(Sj • Tij)^ 



Sj) 



(1) 



where Sj = (sf,sf,sf) are spin vectors of unit length at 
position Ti = (rf,r\,rl) in layer Xi and fij = — fj. 
The positions ?i are normalized such that nearest neigh- 
bors obey ruj\ ~ 1. J is the nearest-neighbor exchange 



coupling constant, D^) and D£' are the local uniax- 
ial anisotropies of second and fourth order, respectively, 
B denotes the external magnetic field with the effec- 
tive magnetic moment [i of the spins incorporated, and 
lu = hq [x 1 1 'Am a 3 is the strength of the long range dipole 
interaction on a lattice with lattice constant a (jxo is the 
magnetic permeability) . 



(4) 
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To calculate the ground-state energy per spin we as- 
sume translational invariance of the spin structure par- 
allel to the film. This assumption is not correct rigor- 
ously since it can be shown that for a perpendicular ori- 
ented magnetization, for instance, a state with striped 
domains is energetically slightly more favorable. How- 
ever the corresponding energy difference for ultra-thin 
films is of order e~ J / 2u and therefore negligible for real- 
istic parameters ptFe- or Ni-films showing spin reorien- 
tation transitional. 

Assuming translational invariance in the xy-plane the 
summation over all spins within a plane can be done ex- 
actly resulting in the energy per surface spin 

J L 

E (s) = £ 

fJ^.U—l 

L / ± 

- 2 £ *V-H V ' \ 
1 n,v=i \ -1 

-£^ 2) (^) 2 + ^ 4) (^) 4 + S-^ (2) 



,sx). The quantities [i and v denote 
„j is the number of nearest neighbors 



with s — (si, . 
layer indices, z 

between layer \m and is, and ^u-i/i are constants arising 
from a partial summation of the dipole interaction- The 
quantities &s have been calculated previouslyOEJ and 
they are listed together with zs in Table |. 

With an external magnetic field B — (0, Bu , B±) in the 
yz-plane, all spins s p are confined to this plane. They 
therefore can be expressed by_ their azimuthal angle 



(0,sini9 p , cos$ M ). Eq. (H) thus can be rewritten as 



E(0) 



1 L 



3o> 



£ [z?f cos 2 ^ + cos 4 ,?, 



By sin $ M + i?^ cos 1?^ 



(3) 



with # = . The ground state is obtained 

by minimizing the energy E($) with respect to In 
zero external field two stationary points of the energy 
given in Eq. (||) are easily identified to be given by i) x = 
(0, . . . ,0) and $>> — (-|, . . . , J), respectively. We define 
a total anisotropy per surface spin in zero field K by the 
corresponding energy difference, K = Effl') — EiJ)^). 
This quantity is given by 



K = 



£ (4 2) 



3u> 
1 



L 



(4) 



The first term is the sum of the anisotropy constants of 
second and fourth order while the second term is due 
to the dipole interaction. Note that this dipole term is 
identical to the dipole anisotropy per unit area -j/iom 2 
calculated within continuum theory, but with additional 
surface correction, as 



3u> 



L 

■ £ ®\»-«\ 

/i,Z/=l 



L 



3u> 



$i+0(* 2 ). (5) 



For K > a perpendicular magnetization is more fa- 
vorable than an in-plane magnetization and vice versa. 
However, in certain parameter intervals additional sta- 
tionary points appear which may lead to an even lower 
energy resulting in a canted spin structure. This will be 
discussed in detail in Section [V| 

In general the minimization of Eq. (|^) has to be done 
numerically. For realistic parameters appearing for in- 
stance for Fe- or Ni films, however, the exchange interac- 
tion is by far the largest term in the Hamiltonian leading 
to a nearly collinear spin structure. In this situation the 
anisotropy terms can be treated as small perturbation 
and as a consequence the minimization can be done to a 
large extend analytically. 



III. PERTURBATION CALCULATION 

We define an averaged angle, 9 = ■r ®v an d devi- 

ations from it, e„, so that = 9 + e v and Ylt=i e ' y = 0- 
Finite e u appear due to the various anisotropy terms and 
they are therefore small for anisotropy terms (including 
the external magnetic field) which are small compared to 
the exchange energy. This will be assumed in the follow- 
ing. Under these circumstances a perturbative treatment 
is possible. We decompose the energy Eq. (||) into two 
parts, 



E0) = E^(9)+SE(9,e) 



(6) 



with SE(9, 0) = and 



E^\e) 



j L 2^ L 



i 



cos 2 ]T £f -cos 4 0]Tl^ 
L{B\\ sm9 + B ± cos 9) 



(7) 



An expansion of SE(9, e) in terms of e then gives 

SE{9,e) = a{9)-e + ^e-C-e + 0{e) 3 (8) 
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lattice 
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Z\ 


ZS>1 


$0 


$i 


$8>1 


sc(OOl) 
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1 





9.0336 


-0.3275 


~ -16^ 2 e- 2 ^ 


fcc(OOl) 
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4 





9.0336 


1.4294 


~ =Fl67r 2 e- v ^' r ' 5 


bcc(OOl) 





4 





5.8675 


2.7126 





Table I: Number of nearest neighbors zg and dipole sums $a for different lattice types. 8 is the distance between layers. 



where we have introduced an obvious matrix notation. 
The gradient 



d 



6E(6,e) 



e = 



is given by 



a(9) = A(0) sin(26>) 



(9) 



(10) 



with 



A x {9) = Df + 2Df cos 2 e-^Y, «W (11) 



Thus, to lowest order the anisotropy terms are linear in 
e while the exchange term expressed in Eq. (g) by the 
matrix C with matrix elements 



(12) 



A=l 



is quadratic in e. 

The minimum of SE(9, e) appears for e v of the order of 
the anisotropy terms showing that the truncated Eq. (||) 
gives the correct energy up to second order in e. Note 
that up to this order the Zeeman term enters only Eq. (Jt]). 
Therefore, at this level of truncation 9 agrees with the 
azimuthal angel of the averaged magnetization. 

It can be easily seen from the definition Eq. ( ^2[ ) that 
e*o = (1, ■ ■ ■ ,1) is an eigenvector of C with eigenvalue 
zero. With this vector it is convenient to rewrite the 
constrain X^=i Ci/ = as a scalar product, e*o • e = 0. 
This notation will be used in the following. 

The minimalization of the energy is done in two steps. 
First we keep 9 fixed and minimize with respect to e„ un- 
der the constraint eb-e = 0. The corresponding energy at 
the minimum, E(6), is accessible for instance by varying 
the external magnetic field and it is precisely this quan- 
tity which for instance is needed to calculate the FMR 
signal. Finally the ground state energy is obtained by 
minimizing E(9) with respect to 8. 

The variation with respect to e v is achieved by intro- 
ducing the function 

y(9,e)=a{9)-e+^e-C-e + \e -e (13) 

where A denotes a Lagrangian multiplier. Stationarity of 
>J(6>, e) gives 



C-e + a(9) + \e a = 6. 



(14) 



Taking the scalar product with eo and noting that e*o-C 
the multiplier A is obtained as 



A = -— e • 3(9). 



Thus e*is determined from 



(15) 
(16) 

with identity matrix 1 and a matrix E with E^ u = 1 
for all matrix elements. To solve this equation for e we 
introduce the pseudo-inverse of the matrix C, which 
in our case fulfills 



L 



E ■ a{9) = 



Ct = 1 - yE. 

Li 



(17) 



The matrix is uniquely defined if one requires that 
it is a symmetric matrix with eigenvector e*o and corre- 
sponding eigenvahp,zero. The matrix elements of are 
explicitly given byLa 



(7t = 1 

2LJ Zl 



L 2 -l 
L + l 



L\ix- 



L + l 



(18) 



It is easy to see that with the help of this matrix Eq. ( |l6| ) 
can be rewritten as 



C • (e + C t • a(9)) = 0. 



(19) 



Since eb is the only eigenvector of C with eigenvalue zero 
the term in brackets has to be parallel to eb- Multiplying 
this term by eb and using eb • e — and eb • = it 
follows 

e = -C^-a(9). (20) 
Inserting into Eq. ( |l3| ) we get the final result 

E{9) = E^(9)+SE{9) (21a) 
5E{9) = -ia(6»)-C t -a(6») + C(e) 3 , (21b) 

where we used the general property = • C ■ of 
the pseudo inverse. The ground state energy is obtained 
by minimizing E(9) with respect to 9. 

Eq. (|H) is the main result of this work, giving a general 
expression for the ground state energy of a thin magnetic 
film in second order perturbation theory. The influence of 
a non-collinear spin structure on the ground state energy 
will be discussed in the following. 



4 



IV. RESULTS 

In the following we drop terms of order 0(e) 3 in E(9) 
and we specialize to a special case in order to obtain 
analytic results. We neglect the exponentially small ef- 
fective dipole interactions between layers with distance 
larger that one, i.e. ^s>i — 0, and we assume that the 
anisotropics D^' which enter the Hamiltonian Eq. (Q) 
are constant within the thin film but may deviate from 
its constant value at the surface (A = 1) and at the in- 
terface to the substrate (A = L), i.e. 



D (n) 

A x (9) 



with 



A s ,m = d. 



(2) 



2Dl 4 ) cos 2 9 



3uj 



It is easy to see that 

1(9) -C^ 1(9) = Ct tl (A 2 s (e) + A?(9)) 
+ 2C\ L A a (9)A i (8) 



(22) 
(23) 



(24) 



(25) 



since -eo = and C\ 1 — C\ L . Then the second order 
correction calculated in the previous section (Eq. ( |21b) )) 
can be written as 



5E(9) = A(0,L)sin 2 (26») 



(26) 



with 



A(9,L) 



L- 1 

" 8Jzi 



±-^(A a (8)+M8)f 



(A s (9) - Ai(6)Y 



(27) 



Note that from now on L can be considered as continuous 
parameter and all quantities are explicitly L-dependent. 
Inserting a(9) and (9) into Eq. © and introducing 
the quantities 

J, 



K (L) = — — (zL — 2z\) 
K 2 (L) = LD& + D P + D 



(2) 



2ttL - -$x ) (28a) 
(28b) 

K 4 (L) = LD^ + D W + Df ] (28c) 
we can finally write for the energy per surface spin 

E(9,L) = K a (L) + A(9,L) sin 2 (29) 

- K 2 (L) cos 2 9 - Ki(L) cos 4 8 

- L(B\\ sin 9 + B ± cos 9) (29) 

Note that the total anisotropy energy K introduced in 
Eq. (|) fulfills 

K(L) = K 2 (L) + K 4 (L), 



as 5E(9) vanishes at the collinear stationary points ?P 
and d 1 - , respectively. K 2 (V) and K^(L) contain the mi- 
croscopic anisotropy parameters and the dipole terms of 
the film averaged over the different layers. 

It is easy to see that an equation for E(9 7 L) in the form 
given by Eq. ( |29| ) often introduced phenomenologicallyu, 
but without the A-term, is obtained if one assumes that 
all spins in the film are strictly parallel. The important 
point to note here, however, is the fact that an addi- 
tional anisotropy energy A(#, L) enters Eq. j2S| ) which is 
connected to non-collinear spin structures originated by 
inhomogeneities in the magnetic film. Indeed, this quan- 
tity only vanishes in the homogeneous case A s = Ai = 0. 
However, for a magnetic thin film the amplitudes A\ 
in general are not constant. Even if the microscopic 
anisotropy constants are homogeneous (which is un- 
likely to occur for a realistic film) this is not the case for 
the dipole term. 

To discuss the implications of this additional 
anisotropy term A(9, L) we first consider the case that 
there is no microscopic uniaxial anisotropy of fourth or- 
der, = 0. In this case A(9,L) — A(L) is indepen- 
dent of 9. Thus for an inhomogeneous distribution of am- 
plitudes A\ , an effective anisotropy term of fourth order 
in cos 9 is generated although there is no corresponding 
anisotropy term of this order in the Hamiltonian. 

If there exists a microscopic anisotropy term of fourth 
order the situation is more complicated: A becomes 9- 
dependent meaning that higher order anisotropy term of 
up to eights order are generated in E(9,L). 

Finally we mention that the quantity A(9 1 L) can be 
further simplified in two common special cases: In the 
case of a symmetric film 
Ai(9) and therefore 



Di n) we get A s (9) 



A i=s (0,L) = - (L ^ 2) Al{9\ (30a) 
while for the case D\ n) = and dP + D ( s 4) »|$ne 



have A s (9) > Ai(9) and 
A i=0 (9,L) = - 



As an important application of these results we will 
study spin reorientation transitions in the next section. 



V. SPIN REORIENTATION TRANSITION 

The direction of the magnetization in the ground state 
for a given thickness L is obtained by minimizing E(9, L) 
(Eq. (§§)). If the total anisotropy energy K(L) (Eq. @) 
changes sign as function of L, a spin reorientation tran- 
sition takes place in which the direction of the magne- 
tization changes either continuously or discontinuously 
depending on the specific form of E(9,L). In the first 
case a so-called canted phase appears. Analytic results 
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for the width and the position of this phase will be de- 
rived in this chapter. 

We decompose K(L) in volume and surface part the 
usual waylij to get 



K(L) = LK V + K s + Ki 



with 



Note that A s s(6) from Eq. fl24|) can be written as 



A sA {6) = K s;i + D^cos26. 



(31) 

(32a) 
(32b) 

(33) 



A spin reorientation transition occurs if the total 
anisotropy energy K(L) passes through zero as function 
of L. If K s + Ki > this happens for sufficiently large 
dipole interaction with increasing L, as then K v < 0. The 
corresponding transition is from perpendicular magneti- 
zation at small L to an in-plane magnetization for large 
L possibly with a canted magnetization in between. This 
type of transition occurs for Fe-films. The opposite sce- 
nario can occur for negative K s + K\ if a positive volume 
anisotropy K v > is present as observed in Ni-films. 
Thus, to lowest order the critical thickness is explicitly 
given by K(L r ) = 0, leading to 



K s + Ki 



For Fe/Ag(100) films + D 

the other quantities £> v 4 \ Dy and $i are negligible 
and we get J* rj 5.5 in good agreement with numerical 
calculations^. 

For L in the vicinity of L r the minimum of E(&) may 
occur at a finite 6, i.e. a canted phase occurs. To deduce 
the limits of stability of the two phases for which 6 = 
and 6 = -| , respectively, we expand Eq. ( p9| ) around these 
angles. From the sign of the corresponding expansion co- 
efficient it follows that in general there are two transitions 
of second order at thicknesses Li and L^~, respectively. 
The phase with 6 = becomes unstable at L'j where 



,(2) 



(34) 



37w. In this case 



K{Ll) + K 4 m) + 4A(0, L|!) = 



(35a) 



at this point. With increasing thickness the parallel 
phase with 6 = ? becomes stable at L^r where 



K(L^) - Ki(L^) - 4A(f ,L±) = 0. 



(35b) 



For K±(L T ) + 4A(L r ) = both transitions coincide re- 
sulting in a jump from 6 = to 6 = ~ at L x . This is 
always the case for L = 1 and in the symmetric case 
also for L = 2 provided vanishes. Otherwise a 

canted phase (K 4 + 4A > 0) or a region with hystere- 
sis (K4 + 4A < 0) appears as described below. Note that 



in the phases 6 = and 6 = f , respectively, e vanishes 
according to Eqs. ( |ic| , po|) showing that in these phases 
all spins are strictly parallel. This is not the case in the 
canted phase. Note also that for finite magnetic fields 
which are neither perpendicular nor parallel to the film 
minimalization of Eq. (|9|) leads to a 6 between zero and 
t[ and therefore to a noncollinear spin structure. 

The difference of the thicknesses at which the two 
collinear phases become instable defines the width 5L T = 
— Lr of the canted region which can be expressed as 



SLr = 2jC 4 (£ r )+4(A(0,L r ) + A(f,£ r )) 
L r K s + Ki 



(36) 



Thus the fourth order anisotropy energy increases 
the width of the canted phase but even without such a 
term a canted region can be observed due to the effective 
anisotropy A(6>, L T ). If the numerator of the right hand 
side of Eq. ( |36|) is positive, a canted phase occur, while 
for negative numerator we find a discontinuous transition 
with hysteresis. 

A similar calculation can be done in finite magnetic 
fields. If the field is orientated perpendicular to the film 
the thickness at which the phase 6 = becomes instable 
is shifted by 



8l} 

L r 



B± 
2K V 



(37a) 



while for fields parallel to the film the corresponding shift 
is given by 



5Lj 



Bn 



Lr 2K V 



(37b) 



A phase diagram for finite temperatures and field has 
been obtained within mean field theory previouslyEj. For 
small external fields the shifts of the phase boundaries 
obtained are linear in the field similar to the present sit- 
uation. 



VI. CONCLUSION 

Starting from a microscopic model the ground state 
energy of a thin ferromagnetic film as function of the 
direction of the magnetization is calculated. Explicit 
expressions for this energy are obtained which contain 
important anisotropy contributions due to non-collinear 
spin structures in certain parameter intervals. The mi- 
croscopic parameters entering the Hamiltonian are not in 
a simple way related to the ground state energy. This is 
important for a comparison of measured and calculated 
anisotropy parameters. Our investigation shows that in 
generally a canted phase is obtained and that the corre- 
sponding transitions into this phase are of second order. 
Analytic expressions are obtained for the width of the 
canted phase and its shift in external magnetic fields. 
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